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Abstract 



We introduce an alternative approach to the third order helicity of a volume preserving 
vector field B, which leads us to a lower bound for the L 2 -energy of B. The proposed 
approach exploits correspondence between the Milnor /2123-invariant for 3-component 
links and the homotopy invariants of maps to configuration spaces, and we provide a 
simple geometric proof of this fact in the case of Borromean links. Based on these 
connections we develop a formulation for the third order helicity of B on invariant un- 
linked domains of B, and provide Arnold's style ergodic interpretation of this invariant 
^ . as an average asymptotic /2123-invariant of orbits of B. 

m 
in 

1 Introduction 

00 

■ A purpose of this paper is to develop a particular formula for the third order helicity on 

certain invariant sets of a volume preserving vector field B. The third order helicity, [3], is 
an invariant of B under the action of volumorphisms isotopic to the identity (denoted here 
by SDiff (M)). Importance of such invariants stems from the basic fact that the evolution 
of the vorticity in the ideal hydrodynamics or of the magnetic field Bq in the ideal magne- 
tohydrodynamics (MHD), occurs along a path t — > g(t) G SDiff (M), P p. 176]. Namely, 
B(t) = g*(t)B which is a direct consequence of Euler's equations: 

|b + i„,b] = o, ±m-v. (i.i) 

One often says that the magnetic field B is frozen in the velocity field v of plasma, and the 
action by SDiffo(M) is frequently referred to as frozen-in-field deformations. A fundamental 
example of such invariant denned for a general class in SVect(M) is the helicity H 12 (-Bi, B 2 ) 
defined for a pair of vector fields B\ and B 2 on M = S 3 or a homology sphere. Helicity has 
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been first introduced by Woltjer, [39], in the context of magnetic fields, and is a measure 
of how orbits of Bi, and B 2 link with each other. This topological interpretation of helicity 
has been made precise by Arnold, who introduced the concept of the average asymptotic 
linking number of a volume preserving vector field B on M, [2]. The subject has been further 
investigated in [H El [25], [36] (see [3] for additional references), where authors approach higher 
helicities via the Massey products under various assumptions about the vector fields or their 
domains (work in [16j E3 [37] concerns yet other approaches to the problem) . Extensions of 
the helicity concept to higher dimensional foliations can be found in [TBI 1331 [23] and recently 
in pig. 

In this paper we present an alternative to these approaches, which is a natural extension 
of the notion of the linking number as a degree of a map, and exploits relations to the 
homotopy theory of certain maps associated to the link. The paper has two parts. 

In the first part we show how the Milnor /i 12 3-invariant: Jx\2z{L) of a parametrized Bor- 
romean link L in S 3 can be obtained as a Hopf degree of an associated map to the configu- 
ration space of three points in S 3 . The presented approach to link homotopy invariants of 
3-component links has been proposed in [20] for n-component links in M 3 , and conveniently 
simplified for 3-component links in S 3 in the joint work [12], where the full correspondence 
between /2 12 3(Z/) and the Pontryagin-Hopf degree is proved. In Section [2] we present the orig- 
inal proof of this correspondence in the Borromean case, which is sufficient for our purposes. 

The second part of the paper is discusses a new definition of the third order helicity 
denoted here by Hi23(-B; T), where B is a volume preserving vector field having an invariant 
unlinked domain T C S 3 . The simplest of such domains are three invariant handlebodies 
in S 3 which have pairwise unlinked cycles in the first homology. This includes the case of 
Borromean flux tubes already investigated in [51 [25] . In Section [5] we develop an ergodic 
formulation of H 12 3 (B;T) as an average asymptotic /2123-invariant, in the spirit of Arnold's 
average asymptotic linking number, which allows us to extend the definition of the invariant 
to topologically more complicated unlinked domains. We also derive, in Section [7J a lower 
bound for the L 2 -energy of B in terms of H 123 T). 
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the valuable input. I am also grateful to my advisor Robert Ghrist who introduced me to 
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link maps and Milnor /2-invariants. The author acknowledges financial support of DARPA, 
#FA9550-08-l-0386. 
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2 The Milnor //^-invariant and the Hopf degree. 

The /1-invariants of n-component links in S 3 have been introduced by Milnor in [28j [29] as 
invariants of links up to link homotopy. Recall that the link homotopy is a deformation of 
a link in S 3 which allows each component to pass through itself but not through a different 
component. Clearly, this is a weaker equivalence than the equivalence of links up to isotopy 
where components are not allowed to pass through themselves at all. The fundamental 
example of a /2-invariant is the linking number (denoted by //12) which is a complete invariant 
of the 2-component links up to link homotopy. In the realm of 3-component links the relevant 
invariants are the pairwise linking numbers p,\2, Ji^z, ^32, and the third invariant /2123 in 
^gcd(/ii2,/i23,Ai32); which is a well defined integer, if and only if, // 12 = //23 = /W32 = 0. In the 
second part of the paper we will interpret this statement as a topological condition on the 
invariant set of a vector field. A precise definition of /i-invariants is algebraic and involves the 
Magnus expansion of the lower central series of the fundamental group: tt\{S 3 — L) of the link 
complement. We refer the interested reader to the works in [28, 29J. In the remaining part 
of this section we will prove that p>i23(L) is a Hopf degree for an appropriate map associated 
to the link L, provided that the link is Borromean, i.e. the pairwise linking numbers are zero 
(note that the Borromean links are more general then Brunnian links, [28]). 

Let us review basic facts about the Hopf degree J$?(f) of a map / : S 3 — > S 2 , (see 
e.g. [6]). A well known property of the Hopf degree is that Jt? : / — ► J4?(f) provides an 
isomorphism between ir 3 (S 2 ) and Z. Recall that up to a constant multiple we may express 
Jjf(f) as (M = S 3 ) 

Jf?(f) = / aAf*v= aAw= / a A da, (2.1) 

J M J M J M 

where v is the area 2-form on S 2 , and a satisfies to = f*v = da. Notice that f*v is always 
exact since the cohomology of S 3 in dimension 2 vanishes. We may also interpret Jff(f) as 
an intersection number, [6j [30] . Namely, consider two regular values p\ and P2 £ S 2 of the 
map /, then Zi = / _1 (pi) and I2 = f' 1 ^) form a link in S 3 , and the integral formula ( 12.1ft 
can be interpreted, as the intersection number of l\ with the Seifert surface spanning I2. 

J?(f) = \k(hJ 2 ) ■ (2.2) 

If we replace S 3 with an arbitrary closed compact orientable 3-dimensional manifold M, we 
may still obtain an invariant of / : M — > S 2 this way, provided that the condition f*u = da 
holds. 

Proposition 2.1. Let M be a closed Riemannian manifold, and v £ Q 2 (S 2 ) the area form on 
S 2 . The formula ( 12.11) provides a homotopy invariant for a map f : M — > S 2 , if the 2-form 
f*v is exact. Up to a constant multiple this invariant can be calculated as an intersection 
number defined in (12.21) . where l\ = f~ 1 {p\) and I2 = f~ 1 {p2) form a link in M, where both 
l\ and I2 are null-homologous. 
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Proof. Given a homotopy F : / x M i-> S 2 , f 1 = F(l, ■ ), f = F(0, ■ ), we define Q = F*v. 
We have 

Cj = F*u = da, LO\ — flv — i\F*u — da\, ujq — f$v — i$ F*v — dato, 

where i : M M x I, io(x) = (x, 0), i\ : M M x J, = (x, 1) are appropriate 

inclusions. Potentials: «|r i xM ) and Q;|| 1 } xM , a i differ by a closed form 

«| {0}xM -ttO = /?0, "|{i}xM _ttl = d(3 = d(3 1 = 0, 

therefore the Stokes Theorem immediately implies that Formula (12.11) is independent of the 
choice of the potential. For the proof of invariance under homotopies we revoke the standard 
argument in 0, p. 228] 

0= / v A v = / CjACj = / uAda 

JF(MxI) JMxI JMxI 



d(u A a) — Cj A a = \ u)\ A ai\ — \ ujq A olq 

Mxl Jd(MxI) JM JM 



da\ A a\ — I dao A«o = J$?(fi) — Jrf?(fo). 
m Jm 

The interpretation of Jt?(f) as the intersection number (12.21) is the same as in [6l p. 230]. □ 

Given a 3-component parametrized link L = {L\, L 2 , L3} in S 3 we wish to associate a certain 
map Fl : S 1 x S 1 x S 1 1— > S 2 to it, and interpret its Hopf degree as the Milnor /2i 2 3-invariant. 
Recall the definition of the configuration space of k points in M: 

Conf fc (M) := {(xi, x 2 , . . . , x k ) G {M) k \ x { ^ x h for i ^ j}. 

As an introduction to the method we review the Gauss formula for the linking number 
of a 2-component link L = {Li,L 2 } in M 3 . Denote parameterizations of components by 
L\ = {x(s)}, L 2 = {y(i}} and consider the map 

F L : S 1 x S 1 Conf 2 (M 3 ) S 2 , L(s,t) = (x(s),y(t)) . 

where r(x,y) = \ZZy\i is a retraction of Conf2(M 3 ) onto S 2 . It yields the classical Gauss 
linking number formula: 

fl 12 (L) = lk(L 1? L 2 ) = deg(F L ), deg(F L ) = / F* L {v), 

Js 1 xS 1 

where v G Q 2 (S 2 ) is the area form on S 2 . Consequently, the linking number lk(Li,L 2 ), also 
known as the Milnor /^-invariant, can be obtained as the homotopy invariant of the map 
Fl associated to L. Observe that homotopy classes [S 1 x S^S* 2 ] are isomorphic to Z and 
deg : F — > deg(F) provides the isomorphism, we also point out that as sets: [S k , Conf„(IR 3 )] = 
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Figure 1: left: yU 12 3 = ±1 and right: /i 12 3 = ±n. 



7r fe (Conf n (.R 3 )), and [S k , Conf n (5 3 )] = 7r fc (Conf n (S 3 )) (c.f. [Til p. 421]). Thus considering 
the based homotopies, in context of the link homotopy of Borromean links, and base point 
free homotopies is equivalent in this setting. 

In [20], [2Tj. authors consider a natural extension of this approach to n-component parame- 
trized links L in M 3 by considering maps F^ : S 1 x . . . x S 1 — > Conf n (M 3 ) and their homotopy 
classes, we refer to this type of maps loosely as link maps, (c.f. [21J). In particular, Kohno J2U] 
proposed specific representatives of cohomology classes of the based loop space of Conf n (IR 3 ) 
as candidates for appropriate link homotopy invariants of L. It has been observed, in [12J, 
that in the 3-component case it is beneficial to consider Conf3(5' 3 ), and L C S 3 , since the 
topology of Conf3(S' 3 ) simplifies dramatically (in comparison to Conf3(M 3 )). We review this 
simplification in the following paragraph as it is essential for the proof of the main theorem 
in this section. 

Consider a 3-component link L = {Li,L 2 ,L 3 } in S 3 parametrized by {x(s),y(t), z(u)} 
and the following map 

F L : S 1 x S 1 x S 1 Conf 3 (5 3 ) S 2 , L(s, t, u) = (x(s),y(t), z(u)), (2.3) 

where we denote by H : S* 3 x M 3 x (R 3 \ {0}) i— > S 2 the projection on the S 2 factor, first 
concluding that Conf 3 (5' 3 ) C S 3 x S 3 x S 3 is diffeomorphic to S 3 x Conf 2 (IR 3 ) = S 3 x M 3 x 
(R 3 \ {0}), and consequently deformation retracts onto S* 3 x S 2 . Considering S* 3 as unit 
quaternions, the map H can be expressed explicitly by the formula, [12J: 



Conf 3 (5 d ) 3 {x,y,z) — >- — - — r — - — r\r^S } 2.4) 

||pr(x 1 ■ y) — pr{x ■ z)\\ 

where • stands for the quaternionic multiplication, _1 is the quaternionic inverse, and pr : 
S 3 — > M. 3 the stereographic projection from 1. As a result one has the following particular 
expression for F L : 

pr(o:(s)- 1 ■ z(u)) - pr(x(g)~ 1 ■ y{t)) 
L{ ' ' ' ||pr(x(s)- 1 -z(M))-pr(x(s)- 1 -y(t))|r K ' 
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At this point we note that one has a freedom in choosing the deformation retraction 
H in (12.41) . but the above particular formula makes the proof of the main theorem of this 
section possible. Let T = S 1 x S 1 x S 1 denote the domain of Fl, notice that, thanks to ( 12. 51) . 
restricting Fl to the subtorus T 23 in the second and third coordinate (t, u) of T, we obtain 
the usual Gauss map of the 2-component link {x~ x ■ L 2 , x~ x ■ L 3 }. Since the diffeomorphism 
x^ 1 - of S 3 is orientation preserving we conclude deg(F L | T23 ) = lk(L 2 , L 3 ). We claim that for 
any 2-component sublink {Li,Lj} of L: 

deg(F L | Tij .) = ±lk(A, Lj), 1 < i < j < 3, (2.6) 

where i,j index the coordinates of T. Indeed, since already true for i = 2 and j = 3, the 
general case follows by applying a permutation a G S3 of coordinate factors in Conf 3 (5' 3 ) C 
(S 3 ) 3 . Notice that a is a diffeomorphism of Conf 3 (S' 3 ) C {S 3 ) 3 either preserving or reversing 
the orientation (which explains the sign in (12.61) ). We infer (I2.6P because a induces an 
isomorphism on homotopy groups of Conf3(S' 3 ). 




Arc A\ is a part of the unit circle on xy-plane, A 2 is a part of circle of radius 1 + e. 



The main theorem of this section is 

Theorem 2.2. Let L = {L x , L 2 , L 3 } be a 3-component Borromean link in S 3 , consider the 
associated map Fl defined in (12. 3p . The Hopf degree of this map satisfies 

Jif(F L ) = ±2fi 123 , (2.7) 

where the sign depends on the choice of orientations of components of L. 

Proof. By the Borromean rings we understand any 3-component link with /2 12 3 = ±1. Every 
such link is link homotopic to the diagram presented on Figure [T] (where the sign can be 
determined from the orientation of components). The proof of the theorem can be reduced 
to the case of Borromean rings L Borr , as follows: if L and V are link-homotopic then Fl 
and Fl> are homotopic maps. By the Milnor classification of 3-component links up to link 
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Figure 3: Four positions corresponding to angles: s = 0, |,7r, ^ small arrows next to iV 
and S indicate a motion as s /. 



homotopy (see |28j), every 3-component link L with zero pairwise linking numbers and 
p,i23 = ±n is represented by the right diagram on Figure [TJ Consequently, up to homotopy, 
the associated map Fl can be obtained from Fi Bori by covering one of the S 1 factors in T, 
n-times. Therefore, in order to prove the claim it suffices to show 

^(F L Borr) = ±2 . (2.8) 

According to Proposition 12. 11 J4?(Fl) is well defined for a link L C S 3 provided F£z/ e Q 2 (T) 
is trivial in H 2 (T), which is true thanks to (12.61) and because the pairwise linking numbers 
of L are zero. The method of proof relies on a direct calculation of Jif(F L Boi), for a carefully 
chosen parametrization of L Bor in S 3 . This calculation is achieved by visualization of the 
link lg t N = U In in T, and application of Formula (12.21) . where 

Is '■= F[ior(S), l N := F~ BoI (N) 

are preimages of the North pole iV = (0,0,1) and South pole S = (0,0,-1) in S 2 C M 3 . 
Notice that [F* L v] = in H 2 (T), if and only if, [l s ] = and [l N ] = in #i(T). 

We begin by identifying S 3 with the set of unit quaternions in M 4 with standard coordi- 
nates: 

(w, x,y,z) = w + x i + y j + z k, 

and choosing a specific parametrization of the Borromean rings L Bor in S 3 . That is, define 
the Li component of L Bor to be the great circle in S 3 through 1 and k, parametrized as 

x(s) = cos(s) + sin(s) k, ( x(s)^ 1 = cos(s) — sin(s) k ) . 

Observe that pr(x(s)) parameterizes the z-axis in R 3 . Figure [2] shows how to define the 
second and the third component {L2,L%} of the Borromean rings L Borr in M. 3 considered as 
an image of S 3 — {1} under the stereographic projection pr : S 3 C M 4 — > M 3 from 1 G S 3 . 
L2 will bound the annuli with a rounded wedge removed, i.e. an arc A\ of the circle of radius 
1. The arc A2 belongs to the circle of radius r e = (1 + e) in the (x, y)-plane. The component 
L3 is chosen to be a vertical ellipse linking with L2. 



7 




Figure 4: Projection of l$ t N on the sw-face and st-face of T. The strands Is (solid line) and 
In (dashed line) are oppositely oriented since Is and In are null homologous in T. 

Next we focus on the Formula ( 12. 5ft . observe that multiplication by x(s) _1 has an effect of 
a rotation by angle s in (w, 2;)-plane and (x, y)-p\ane of M 4 , which can be directly calculated: 

x(s)^ 1 ■ (w, x, y, z) = (cos(s)n; + sin(s)z, cos(s)a; + sm(s)y, 

cos(s)y — sm(s)x, cos(s)z — sin(s)u>) . 

The flow defined by this ^-action is tangent to the great circles of S 3 , thus the projected 
flow on M. 3 , via the stereographic projection pr, presents the standard picture of the Hopf 
fibration. Let us call an invariant Hopf torus an r-torus, if and only if, it contains a circle 
of radius r in the (x, y)-plane. Without loss of generality we assume that L2 on Figure [2] 
belongs to the r e / 2 -Hopf torus. Every point on a r-torus traces a (1, l)-curve under the S 1 - 
action. For sufficiently small e, this motion can be regarded as a composition of the rotation 
by angle s in both the direction of the meridian and the longitude of a r-torus. Therefore, 
for different values of s the ^-action "rotates" the components L 2 and L 3 , by sliding along 
the Hopf tori by angle s in the meridian and the longitudinal direction. We denote resulting 
link components by 

L s 2 (t) = pr(x _1 (s) • y(t)), and L s 3 (u) = pr(ar _1 (s) • z{u)). 

The unit circle on (x, y)-plane is left invariant under this action and therefore can be consid- 
ered as the "axis of the rotation" . This justifies the choice of the particular shape of L Borr 
pictured on Figure[2j Next, we seek to visualize the projection of Is,n = IsUIn on the sw-face 
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and st-face of the domain T of F l b oi parameterized by (s, t, u) G T, (it is convenient to think 
about T as a cube in (s, t, u)-coordinates, see Figure @J. For example when s = 0, (x(0) = 1), 
a point (0,to,Uo) belongs to In, if and only if, the vector v = pr(y(t )) — w(z(uo)) points 
in the direction of iV = (0, 0, 1), analogous condition holds for direction S = (0, 0, —1), and 
Is- In order to determine a diagram of Is,n, we must keep track of the "head" and "tail" of 
the vector v s = L s 2 (t) — L^(u), for various values of s and record values of t and u for which 
m s points "North" and "South" (Figure [3]). This reads as the following condition 

(s,t,u) G Is,n, if an d only if, v s || S or N. (2.9) 

Without loss of generality we assume that L s 2 is parametrized by the unit t-interval, and L| 
is parametrized by the unit it-interval. The process of recording values of u and t such that 
(12.91) holds is self-explanatory and is shown on Figure [2] for values s — 0, |,7r, which is 
sufficient to draw projections of 1$ and In on st- and tw-faces of T. Collecting the information 
on Figure [3], we draw the projection of Is,n on the sw-face of T represented by square (A) in 
Figure HI Analogously, the projection of Is,n on the st-face of T is obtained and pictured in 
square (B). In order to obtain the diagram of Is,n we resolve the double points of Diagram 
(A) into crossings. For example, let us resolve the "circled" double point on (A), which 
occurs at s — |- in the left two stands of Is,n- It suffices to determine the value of the 
t-coordinate at this point. Diagram (B) tells us that Is is below In, because T is oriented so 
that the t-axis points above the sw-face (see Figure H]). Resolving the remaining crossings in 
a similar fashion leads to a diagrams of 1$,n presented in squares (C) and (D). Clearly, the 
linking number of Is and In is equal to ±2 in Diagram (C), (as the intersection number of e.g. 
Is with the obvious annulus on Diagram (C)). This justifies A2.8I) . and ends the proof. □ 

Results of Theorem 12.21 and Proposition 12 . II combined with Formula (12. ip allow us to express 
fli23{L) of a Borromean link L as 



//i 23 (I) = / F[uAa= / L*lu A a, 
J J" 



(2.10) 

for Flu = da, u = H*v G fi 2 (Conf 3 (^ 3 )), 

where v is the area form on S 2 , and H : Conf 3 (S' 3 ) — ► 5* 2 is the deformation retraction (as 
e.g. in (I2.4p ). Alternatively, we may view u as a 2-form on (S* 3 ) 3 which is singular along 
the diagonals A C {S 3 ) 3 , and the singularity is of order 0(r 2 ), where r is a distance to A. 
Consequently, u is integrable but not square integrable on (S 3 ) 3 . 

Remark 2.3. Notice that the integral formula (12.101) exhibits the following property of p-123' 
#123 (Li,L 2 ,L 3 ) = sign(o-) / u 123 (^ (T (i),^<7(2),^<7(3)) ) G S 3 . 
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3 Invariants of volume preserving flows. Helicities. 



Given finitely many volume preserving vector fields B±, B 2 , . . . B k G SVect(M) on M = S 3 
or a homology 3-sphere one seeks quantities \(B\, B2, ■ ■ ■ , Bk) invariant under the action of 
volumorphisms isotopic to the identity g G SDiffo(M), commonly known as helicities or 
higher helicities: 

\{B 1 ,B 2 ,...,B k ) = \{g*B 1 ,g*B 2 ,...,g st B k ), for all g G SDiff (M), (3.1) 

where g* is a push-forward by a diffeomorphism g. To distinguish the case of a single vector 
field B (i.e. B — B\ — . . . — B k ) we often refer to \{B) = \(B, B, . . . ,B) as self helicity. 
We elucidated in the introduction a fundamental example of such invariant is the ordinary 
helicity H(Bi,B 2 ) of a pair of vector fields. In the remaining part of this section we review 
well known formulations of the helicity, which will later help us to point out analogies to 
the proposed formulation of the 3rd order helicity. Let T = T\ U T 2 represent two invariant 
sub domains (not necessarily disjoint) under flows of Bi and B 2 in S 3 and let 

T = T 1 xT 2 C Conf 2 (M) C M x M. 

Recall that the formula for H(Bi, B 2 ), from [T9|. |38j, specialized to invariant subdomains 
T = T\ U T 2 may be expressed as 



H 12 (B U B 2 ) = to A L Bl fJ. A L B2 fi, (3.2) 



'TixT 2 

where uo is known as the linking form on M x M. When T = MxM this formula is equivalent 
to a more commonly known expression: H(Bi,B 2 ) = § m lb 1 [i A d _1 (iB 2 A i ) (because u also 
represents the integral kernel of d~ l ). Philosophically, H(Bi,B 2 ) can be derived (c.f. [2]) 
from the linking number of a pair of closed curves, which is expressed by Arnold's Helicity 
Theorem. For orbits {^ >1 (x), 2 {y)} of B\ and B 2 through x, y G M, we introduce the 
following notation for the long pieces of closed up orbits 

= I < t < T} C ^(x), i = 1,2 

^(x) := ^(x) U a(x, $*(x, T)). 

where a(x,y) denotes a short path, [38], connecting x and y in M (see Section Ej). Para- 
phrasing [2 J we state (for the proof also see [38J ) , 

Theorem 3.1 (Arnold's Helicity Theorem, [2]). Given B\,B 2 G SVect(M), The following 
limit exists almost everywhere on M x M : 

fh BlB2 (x, y) = lim ^- lk(^{x) x & 2 T {y)) . (3.4) 

Moreover, m BxB% is in L l {M x Af), and 

Hi 2 (B 1 ,B 2 )= rh Bl B 2 (x,y) /i(x) A fi(y) . (3.5) 

J MxM 
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The function ffiB l B 2 represents an asymptotic linking number of orbits {^ >1 (x), €? 2 (y)}, and 
the identity (I3.5P tells us that the helicity H 12 (B 1 ,B 2 ) is equal to the average asymptotic 
linking number. In coming paragraphs, we will demonstrate, how this philosophy is applied 
to obtain the asymptotic /ii23-invariant for 3-component links and the third order helicity. 

4 Definition of "/2123-helicity" on invariant unlinked han- 
dlebodies. 

In this section we apply the formulation of the /ii23-invariant for the 3-component links in 
S 3 , obtained in Section [2, to define the third order helicity of a volume preserving vector 
field B on certain invariant sets T of B in S 3 . In the following paragraphs as a "warm-up" 
to a more general case treated in Section [5j we consider the case of three disjoint unlinked 
handlebodies: T = TiUT 2 U% in S 3 each of genus g(7i). Henceforth, we use "unlinked" to 
mean "with pairwise unlinked connected components". When T represents three unlinked 
tubes (also known as flux tubes [3]) in M 3 , the third order helicity has been developed by 
several authors [5j [2TJ, [25] via Massey product formula for the /2123-invariant, we compare 
our approach to these known works in Section [HJ 




Figure 5: The simplest unlinked handlebodies: flux tubes T modeled on the Borromean 
rings (left), and unlinked genus 2 handlebodies (right). 

Assume %s have smooth boundary and B to be tangent to d%, we set 

Bi'=B\ z = l,2,3, 

and denote the flow of B on S 3 by $, and flows of restrictions Bi by <3>\ Clearly, such T is 
an invariant set of B. Given any domain T with three connected components {%} we may 
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always associate a product domain in Conf 3 (S' 3 ) as follows 



T := T x x T 2 x T 3 c Conf 3 (S 3 ) c 5 3 x 5 3 xS 3 . (4.1) 

Notice that T is a domain with corners in Conf 3 (5' 3 ), and we use the same notation for the 
product of % as for the union in S 3 . Wherever needed, we also assume that (S 3 ) 3 is equipped 
with a product Riemannian metric. Let a domain T defined in (14. II) . where %r\Tj =0, i ^ j 
and each % is a handlebody in S 3 be called unlinked handlebody, if and only if, the 2-form 
uj G f2 2 (Conf 3 (S' 3 )) defined in Equation (12.101) is exact on T, i.e. uj has a local potential 

uj = da w . (4.2) 
Denote a volume preserving vector field B and an unlinked handlebody T as a pair (B; T). 

Remark 4.1. Since uj is a dual cohomology class to the S 2 factor in Conf 3 (5' 3 ) = S 3 x S 2 , 
uj does not admit a global potential. 

Because each handlebody % has a homotopy type of a bouquet of circles, there is a natural 
choice of the basis for Hi{T^) which consists of cycles {L\ }k=i,..., g (dTi) corresponding to the 
circles. We have the following practical characterization of unlinked handlebodies: 

Lemma 4.2. T is an unlinked handlebody, if and only if, 

1^,^ = 0, for all i, j, i^j (4.3) 

Proof. The standard integral pairing, [6], H 2 (T) x H2CT) — > R implies that a closed /c-form 
is exact, if and only if, it evaluates to zero on all /c-cycles of the domain. By the Kiinneth 
formula H 2 (T) = H 2 {T l x T 2 x T 3 ) is generated by L\ ® U p and by (1231) : 

\\{Llm=\ uj i^j. 



L*fy.L r A 



Therefore the condition (14.31) is necessary and sufficient for uj to be exact on T. □ 
We define the /2i2 3 -helicity of [B\ T) denoted by Hi 23 (-E>; T) or Hi 23 (-Bi, B 2 , B 3 ) as follows: 



H 123 (B; T) = H 123 (B 1 , B 2 , B 3 ) d = J (o^ A da^) A l Bi ^\ A i B ^2 A t B3 // 3 , 



(4.4) 



where [ii denotes the pull-back of the volume form fi on S* 3 under the projection 

rti : S 3 x S 3 x S 3 — > S 3 , ix i (xi ) x 2) x 3 ) = x h (4.5) 

and is a contraction by a vector field Bi. Our notational convention is to denote by ib^i 
both the forms on the base of i\i and the pullbacks: 7r*(i B .yUj). Notice that \i = yUi A \i 2 A /i 3 
is a volume form on the product: S 3 x S 3 x S 3 . There are obvious analogies between 
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Formula (14. 4p above, Formula (13.21) for Hi2(-B;T), and the integral formula f 1 2 . 1 j) for the 
/2 123 -invariant. The 3-form: 

Iuj ■= ol w A da w = Alo, 

plays a role of the linking form as u in Formula ( 13.21) . The main motivation behind definition 
( I4.4p is the ergodic interpretation of Hi2 3 (i?;T) as an average asymptotic /ii2 3 -invariant of 
orbits of B, which will become apparent in Section [51 Formula ( 14.41) can be also regarded as 
the third order helicity of three distinct vector fields Bi, supported on the handlebodies %. 
In Section [7J we indicate how to construct the potential a w from the basic elliptic theory of 
differential forms. 

Theorem 4.3 (Helicity Invariance Theorem). On every unlinked invariant handlebody T in 
S 3 , Hi 23 (-B; T) is 

(i) independent of a choice of the potential a u , 

(ii) invariant under the action of SDiff (S' 3 ), i.e. for every g £ SDiff (S' 3 ): 

H 123 {B ] T) = H 123 {g m B;g{T)) . (4.6) 

Proof. To prove (i) observe for every other potential a'^ of u, the difference (3 = a u — a' w is 
a closed 1-form on T (since u = dot^ = da'J). Therefore, 

H 123 (B;T) -H' 123 (B;T) = J (0 A u) A t Bl ^i A ^2 A t B3 ^ 3 

d{f} A A L Bl yL X A i B2 [i 2 A t Bs fi 3 ) 

( 2 ) 

(3 A a w A i Bl [i\ A L B2 fi 2 A is 3 /i3 = 0, 

JdT 

where in (1) we applied d(i Bi [ii) = (since B^s are divergence free), and in (2): 

tBiH-i\ dT = 0, 

(because each vector field Bi is tangent to the boundary d%), where 

dT = [&T X x T 2 x T 3 ) U (Ti x <9T 2 x T 3 ) U (T x x T 2 x <9T 3 ). 

The proof of (n) is in the style of [27], but adapted to our setting. For any given g £ 
SDiff(S' 3 ), by definition, there exists a path t — ► g(t) £ SDiff (5' 3 ), such that 

g(0) = id 5 3, g(l) = g . 

Denote by V the divergence free vector field on S 3 , given by V(x) = 4:g(t,x)\ t=0 , i.e. g(t) is 
a flow of V, and push-forward fields Bi by 

B\ := ^t)*^ . 
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It is well known (see Appendix |Aj or [15, p. 224]) that 2-forms: L B tfi are frozen in the flow 
of V, i.e. 

^(g(t)*t Btf i) = (d t + C v y B jLi = 0. (4.7) 

We also have a path g(t) = (g(t), g(t), g(t)) in SDiff (S' 3 x5 3 x S 3 ), which analogously leads 
to the vector field V = (V, V, V). (Recall that a tangent bundle T(S 3 ) 3 has a natural product 
structure). Equation (14.71) implies 

(d t + C v ) = (d t + C^LBtfM = . (4.8) 

(In the second equation we merely revoke our notational conventions: LB*m = ti* U g ^* B .fi)). 
Let T(t) = g(t)(T(0)) C Conf 3 (S 3 ), we must show ±H 123 (B[, B\, B\) = 0. Notice that for 
small enough e and t G (t — e, to + e ) we can assume, by (i), that a w is a time independent 
potential obtained from slightly bigger domain T which deformation retracts on T(to), and 
satisfies 

T(t) C T, for t E (to-e,t Q + e) . 
Without loss of generality set t = 0, and ^(0) = id( S 3)3, at t we calculate: 

d d C 

■^\Hi2 3 (g(t)*B 1 ,g(t)*B 2 ,g(t)*B 3 yj = — J a u A da^ A i B t\i x A A B */i 2 A L B t/j, 3 

^-Mt)* ((Xu A dau A L B t Ui A L B t U 2 A L B t LL 3 ) 
T(0) at ' 

A (faj) A L B tUx A t B t/i 2 A t fl t/x 3 , 

T(0) 

where in the last identity we applied (14.81) and the product rule for the Lie derivative. Now 
because u A a u is time independent (for t <E (t — e,t + e)), Cartan magic formula yields 

Zft+viatu) A daj) = C v (a w A da^) = t v d(a w A da w ) + d(L v (a UJ A dajj) 
= d(c v (a UJ A daw)), 

where d(a u A da u ) =wAu = 0. Since B\ are tangent to the boundary of %(t), the same 
argument as in the proof of (i) shows that the right hand side of the previous equation 
vanishes. □ 

Remark 4.4. Notice that the above argument indicates that if we replace A da w by 
virtually any closed 3-form r\ on Conf 3 (S' 3 ) (or (5 13 ) 3 ) we obtain some invariant under frozen- 
in-field deformations. If rj is exact we obtain trivial invariants, therefore the only sensible 
candidates here are cohomology classes of Conf 3 (S' 3 ) = S 3 x S 2 . In dimension 3 it leaves 
us with a dual to the S 3 factor in Conf 3 (5' 3 ). Based on the considerations in Section [2] one 
may argue that an invariant obtained this way is trivial. Indeed, the cohomology class rj 
evaluated on any 3-torus obtained from a 3-component link in S 3 via the map L in (I2.3P is 
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zero. Therefore, one could apply the ergodic approach of Section [5] to show that ft) A Lb^x A 
i^b 1^2 A lbP-3 defines a trivial invariant. The crucial obstacle in extending the formula in (14.41) 
to encompass the whole (S 3 ) 3 is the fact that the potential a u cannot be globally defined on 



The following statement is often seen in literature [SlEj: 

Helicity measures the extent to which vector fields twist and coil around each 



A beauty of Arnold's ergodic approach to the helicity Hi2(-B) is that it makes this statement 
precise, by interpreting H 12 (-B) as an average asymptotic linking number of orbits of B. But, 
it also has a practical application as it allows us to extend our approach to certain invariant 
sets of B. In this section we apply this philosophy to our newly defined invariant W\2z{B; T), 
and interpret it as the average asymptotic /2123-invariant of orbits of B in T . Moreover, this 
ergodic interpretation leads us to an alternative, more intuitive proof of Helicity Invariance 
Theorem 14.31 

We begin by observing that given a volume preserving vector field B on. M and its flow 
$f , we may regard B as three vector fields on (M) 3 . Thus, $, $) induces a natural M 3 
action defined as follows: 



Observe that $ is a volume preserving action on (M) 3 . Our analysis is rooted in techniques 
developed in [21 [Ml [23 EH] , the main tool is the following 

Theorem 5.1 (Multi-parameter Ergodic Theorem, [1]). For any real valued L l -function F 
on (M) 3 ; the time averages under the action in (\5.l\i : 



Conf 3 (S 3 ). 



5 The ergodic interpretation of Hi2s(B]T) 



other. 




T 



I 



r 



I 



■T 



F(x, y, z) = lim 



F($(x, s), $(y, t), u)) ds dt du 



converge almost everywhere. In addition, the limit junction F satisfies 



(V II-FIIl^a/) 3 ) < ||F|| L i(( M )3) ; 
(ii) F is invariant under the action, 
(Hi) if (M) 3 is of finite volume then 




(5.2) 
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Definition 5.2. Define invariant unlinked domain T of B as an arbitrary ^-invariant set, 
with topological closure T which belongs to a larger product of open sets T = T\ xT 2 x %, in 
Coni^S 13 ), satisfying the following 

(A) T admits a short path system S, 

(B) Equation (14. 2ft holds on T , 

where by a system of short paths on T, f3R \B$, we understand a collection of curves S = 
{ai(x,y)} on each open set % such that 

(a) for every pair of points x, y G % there is a connecting curve o~i(x,y) : I i— > % in S, 
0i(O) = x and crj(l) = y, 

(b ) the lengths of paths in S are uniformly bounded above by a common constant. 

Topologically, every ^-invariant set is a union of products of orbits of B in (S" 3 ) 3 . It is often 
convenient to think of the orbits ^-action as a foliation of {S 3 ) 3 . Then ^-invariant sets 
are just union of leaves of this foliation. A fundamental example of an invariant unlinked 
domain is the case of ^-invariant set T contained in the product T = T1XT2XT3 of disjoint 
open unlinked handlebodies %. Note that in this case we do not require B to be tangent to 
d%, and T always admits a short path system as we describe in the following 

Remark 5.3. In [38], Vogel shows that on a closed manifold M, geodesies always provide 
a short path system. When T is contained in the product of unlinked handlebodies T we 
may easily construct such system on T as follows Because % are proper subsets of S 3 we 
generally do not want to use ambient geodesies from S 3 as they may not lie entirely in %. 
To obtain S one puts an artificial Riemannian metric on each % which makes d% totally 
geodesic, and choose S to be geodesies on such Riemannian manifold. Observe that applying 
a diffeomorphism g G Diff(T) to S results in the system gS on g(T). 

The following result is an analog of Arnold's Helicity Theorem in our setting, 

Theorem 5.4 (Ergodic interpretation of H 12 3(-B; T)). Given (B;T), the following limit 
(asymptotic [1123-invariant of orbits) exists for almost all (x,y,z) G T: 

m B (x,V,z)=Vim ±Ji 123 (0*(x), 0*(y), 0* 3 (y)). (5.3) 

Moreover, 

H 123 (B;T) = J m B (x,y,z)(i(x) Afi,(y) Afi{z) . (5.4) 
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Proof. The proof is similar to the one in e.g. [38] • Before we start, we must point out the 
following identity (valid for any 3-form (3 on M x M x M and vector fields Bi,B 2 , B% on M) 

(L B3 LB 2 tB x (3) A /il A fl 2 A fJ>3 = /3(-Bl, £ 2 , .63) /Xi A /i 2 A /i 3 ^ 

= /3 A A A 6^3/13 • 

The first equation follows from the definition, the second one is a consequence of the fact 
that lb is an antiderivation i.e. 

L B (a A 0) = (l b o) A f3 + (-l) |o| a A (i B 0), (5.6) 

and ^/ij = 0, for i ^ j (see Appendix |A|) . As a result, 

H123OB1, #2, B 3 ) = J (a. A da^) A A Ib 2 // 2 A Ib 3 // 3 

(i Ba i B2 L Bl (a u A da w )) /ii A ii 2 A /i 3 
m Bu B 2 ,B 3 (x,y,z) iii A /i 2 A/i 3 , 



r 



r 



where m BljB2tB3 := A da w (Bx, B 2 , B3). For convenience, let us set (see (13. 3ft ): 

#(x,y,z;T) := 0*{x) x x tf*(y) . 



Observe that if the orbit ff(x,y,z;T) is nondegenerate, it represents a Borromean link and 
we may apply Formula (12.101) to get 

fii2 3 (^T 1 ^)^T 2 (y)^T 3 (y)) = [ «,Ad« u 

Jff(x,y,z;T) 

a u A da^ + (I), 

[x,y,z;T) 

where the term (I) involves integrals over short paths in S, (see Appendix [B]). For degenerate 
orbits (such as fixed points etc.), the above formula still makes sense because /2i 2 3 is a 
homotopy invariant of the associated map as proven in Section [21 For every (x, y, z) : 

D$ {x ^ z) [di] = B u 1 = 1,2,3, 

where d\ := d s ,d 2 '■— dt,ds := d u , thus we obtain 

*(*,»,*) A daj) = (oiu A da^B-L, B 2 , B 3 )(^ 1 (x, s), $ 2 (y, £), $ 3 (,2, u)) dsAdtA du 
= ^Bi,b 2 ,b 3 (^ 1 (^ ) s), $ 2 (y,t), $ 3 (^, m)) ds Adt Adu . 
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Therefore, 

/ a UJ Ada u) = / / / rn Bl)B2iB3 (^ 1 (x, s),$ 2 (y,t),$ 3 (z,u)) ds A dt A du . 

J ff(x,y,z;T) JO JO JO 

Function m BljB2jB3 is smooth bounded on T and hence L l . Because short paths do not 
contribute to the time average (see Appendix IB1 : 

lim -L(/) = 0. 

Theorem 15.11 applied to the function m Bl)B2jB3 yields almost everywhere existence of the limit 
(I5.3p . Hence we obtain the invariant L 1 -function rh B := fh BljB2tBi on T. The identity (15.41) 
follows from (Hi) of Theorem 15.11 □ 



Theorem 5.5 (Helicity Invariance Theorem-ergodic version). On every unlinked domain T 
in S 3 , (i) and (ii) of Theorem \4-3\ hold for Hi2s(B;T) . 

Proof. The proof of (i) immediately follows from independence of the limit (15.31) of the choice 
of the potential a w . For the proof of (ii) we must show the following 

/ ^Si,B 2 ,-B 3 fJ'l A 1^2 A /i3 = / fix , 9H ,B 2 ,cf*B3 A*i A p-i A 1^3 ■ 

JT Jg(T) 

Theorem 15.11 tells us that m BljB2tBs and rn gitBltgtrB2 ^ tBz admit Z^-averages rh BljB2tB;i and 
^ l g f B 1 ,g t B 2 ,g,B 3 , under actions of Bi and g*Bi respectively. It suffices to show the following 
identity 

fh BljB2jB3 (x,y,z) = rhg. Bl ,g.B 2 ,g.B 3 (g(x),g(y),g(z)), a.e. (5.7) 

then (14.61) is an immediate consequence of Equation (I5.2p . change of variables for integrals, 
and the fact that g preserves volume (i.e. /ij = g*[ii)- Borrowing notation from the previous 
theorem set 

g0(x,y,z;T) := g(0*(x)) x g(^(y)) x g(^(z)) 

= ^ Bl (g(x)) x ff^(g(y)) x 9 T * B \g{z)\ 

where the second identity is a consequence of the fact that the flow of g*Bi is obtained as 
a composition of g and the flow of Bi. Since g is isotopic to the identity, 0(x,y, z;T) and 
g0(x, y, z; T) are homotopic as link maps (for nondegenerate orbits they are in fact isotopic 
Borromean links in S 3 ) and by theorems of Section [21 we have 

AmC^fo V, z; T)) = p, l2 -i(gd(x, y, z; T)). 

As a result of the above identity and (15. 3p we derive a.e. 

rh BltB2tB3 (x,y,z) = lim —p 123 (^(x,y, z;T)) 

= lim —fi 123 (g0(x,y,z;T)) = m g , Bug » B2tg « B3 (g(x),g(y),g(z)). 

I — »oo 1 
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Notice that in the last equation we used the "pushed forward" short paths system: gS. Since 
the lengths of paths in gS are bounded as well, they do not contribute to the limit. This 
proves the identity ( 15. 71) . and consequently ( 14.61) . □ 

Notice that the above argument does not require the Stokes Theorem, and as such may lead 
to further generalizations. Clearly, for H123 (B;T) to be nontrivial T must be of nonzero 
measure. 



6 Flux formula for Hi 2 s(B]T). 

The following formula is a well known property of the ordinary helicity Hi 2 (B; T) of the flux 
tubes T modeled on a 2-component link L = {L 1; L 2 } (see e.g. [2U E]) 

H 12 (5; T) = H 12 {B U B 2 ) = Ik(L ls L 2 ) F\mx(B 1 ) Flux(5 2 ). (6.1) 

Here we show an analogous property for Hi 23 (B; T), when T is an invariant unlinked han- 
dlebody. Recall that {£f }fc=i,..., 9 (,9T;) denotes the basis of Hi(T) defined in Lemma l4~2l 

Proposition 6.1. H 123 (Bi, B 2 , B 3 ) on invariant unlinked handlebodies T satisfies the fol- 
lowing formula 

H l23 (B u B 2 ,B 3 ) = J2fli2 3 (L[,LiL k 3 ) Flvx^W Flux^.{B 2 ) Flux^ k {B 3 ), (6.2) 

where {L\ <8> L 2 <g> L 3 } is a basis of H 3 (T), Flux(Bi) stands for the flux of Bi through a cross 
sectional surface of%, which represents the homology Poincare dual of Lf in H 2 (%, 8%). 

Proof. Recall that the flux Fluxs fc (-B«) of a vector field Bi though a cross-sectional surface 
S fc in % is given by: 

Flux Sfe (5i) = / lb^= \ hkAiBi/J. 



(6.3) 

/ t Bi h k An= / h k (Bi), 



where 1-forms h k represent cohomology Poincare duals of and we applied (15.61) in the 
third equation. For every closed curve 7 C T h k satisfies 

h k = #(7, = deg(7, L k ), 

where deg(7, L° k ) measures how many times 7 "wraps around" the cycle L 3 k . 

For simplicity, we first assume that T is modeled on a Borromean link L = {L\, L 2 , L 3 } 
(such as T Borr on Figure [5]). Then, h = hi A h 2 A h 3 G il 3 (T) is a cohomology class dual to 
the cycle Li ® L 2 <g> L3 in H 3 (T). Define 

_£/~ := L B3 iB 2 i B 1 h = h(Bi, B 2 , B 3 ), 
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which is a smooth function, and let H be the time average of H as in Theorem 15. II It suffices 
to show 

fh Bl ,B 2 ,B 3 = P123 (£) H, a.e. (6.4) 

Analogously, as in the proof of Theorem 14.31 Equation (\QA\i immediately implies Formula 
(16. 2p . Assume the notation of Theorem 14. 3[ for given T consider 0(x,y, z;T). Thanks to 
Theorem 12.21 (see the first paragraph of the proof) we have 

p 123 {0(x,y,z;T)) = # 123 (L) (deg(^(x), Ly) deg(^(y), L 2 ) deg(0*(x), L 3 )) 

= #12 3 (£) / ^ • 



Therefore, 



rn Bl ,B 2 ,B 3 {x,y,z) = lim —/i 123 (^(x,y,z;T)) 

1 — >oo J 

= # 123 (L) hm — / H 



{x,y,z;T) 



/j ll23 (L)H(x,y,z), 



where the last equality is again the consequence of short paths not contributing to the limit. 
From the product structure of T and (Hi) of Theorem 15.11 we get 

H= [ H=(l h l {B 1 )){f h 2 (B 2 )){[ h 3 (B 3 )), 

which combined with Equation (16.31) for fluxes concludes the proof in the case of Borromean 
flux tubes. The proof in the case of a general handlebody is analogous, once we show the 
following 

Lemma 6.2. Let O = {0\, 2 , O3} be a 3-component link in S 3 such that Oi C % for each 
i. Then O is Borromean and 

fii2 3 (0) = h2s(L{, Ll L3) deg(O u L\) deg(0 2 , V 2 ) deg(0 3 , L\) . (6.5) 

i,j,k 

Proof. Thanks to the interpretation of the /2123-invariant in Section [2], it is not only a link 
homotopy invariant, but also a homotopy invariant of the associated map Fo defined in (12. 3p . 
Observe that each component Oi can be homotoped inside of its handlebody % to become a 
bouquet of circles Oi = S 1 V S 1 V . . . V S 1 so that each factor in Oi is a multiple of the cycle 
represented by {L]}in H\(Ti). As a result, we obtain the associated map Fq, and 

Interpreting J4?(Fq) as the intersection number and summing up intersection numbers we 
conclude (16.51) . □ 
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□ 



In the case of Borromean flux tubes, Formula (16. ip reduces to 

H 123 (5i,B2,£ 3 ) =/ii 2 3(^)Flux 2l (S 1 )Flux S2 ( J B 2 )Flux E3 ( J B 3 ), (6.6) 

where Ej denotes homology Poincare duals to Lj in H 2 (%,d%). 

Since the fluxes are invariant under frozen-in-field deformations, Formula ( 16.1 ft is yet 
another proof of Theorem 14.31 in the setting of invariant unlinked handlebodies. In [25] 
the authors develop the same formula for the Borromean flux tubes. This clearly must be 
the case, as we work with the same topological invariants of links via a different approach. 
Additional advantage of our formulation is that we do not have to separately deal with null 
points of vector fields as in [25] . 

7 Energy bound. 

In this section we indicate how the quantity H 123 (B;T) invariant under frozen-in-field de- 
formations provides a lower bound for the L 2 -energy E 2 (B) of a volume preserving field B 
on M = S 3 . We restrict our considerations to the case of an invariant unlinked handlebody 
T, defined in Section HI For the notation used in this section see Appendix 
Recall the definition 

E 2 (B)= [ \B\ 2 = ||S||| 2(M) . (7.1) 

J M 

The ordinary helicity Hi2(-B) provides a well known lower bound (see [HI p. 123]): 

Uh 12 (b)\<e 2 {b), 

where Ai is the first eigenvalue of the elliptic self adjoint operator *d : Q 1 (M) 1— > Q 1 (M) 
(known as the curl operator), * denotes the Hodge star. Importance of such lower energy 
bounds, stems from an area of interest in the ideal magnetohydrodynamics, [32J, as this 
constrains the phenomenon of "magnetic relaxation" , [2] . A need for higher helicities can 
be justified by the fact that one may easily produce examples of vector fields B for which 
H 12 (-B) vanishes, but the energy of the field B still cannot be relaxed. For example, consider 
a classical case of Borromean flux tubes T Borr with B smooth and vanishing outside the 
tubes. Furthermore, assume that orbits of B are just "parallel circles" inside each tube. By 
taking B = Bi + B 2 + B 3 , bilinearity of H 12 ( • , • ) on SVect(S' 3 ) and disjoint supports of B^ 
yield 

3 

H 12 (B) = H 12 (B, B) = + Hl2 ( 5 - B i) = °' 

i=i i^j 

where "cross-helicities" H 12 (-Bj, Bj), i ^ j vanish by Formula (I6.ip . and self helicities H 12 (-Bj) 
vanish because the average linking number of orbits is zero (orbits are just parallel circles). 
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Nevertheless, Formula (16. 2p tells us 

H 123 (£;T Borr )^0, 

as a result we may regard Hi23(-B; T) as a possible "higher obstruction" to the energy relax- 
ation or the third order cross-helicity of B on T . 

To obtain a lower bound for E 2 (B) in such situations we notice that H 123 (S;T) is the 
L 2 -inner product of the 6-forms: *(wAa u ) and Lb^x A £b 2 // 2 A £b 3 jU3- Indeed, from (14.41) we 
have 

H 123 (.B; T) = (*(a w A u), Lb^x A Ib^i A fc^a) L 2 Q 6 ({5 3 )3) • 

Let CV = ||ot w Ao;||l2, for a fixed Riemannian product metric on (S 3 ) 3 this constant depends 
on the domain T in (S 3 ) 3 . We estimate using the Cauchy-Schwarz inequality 

|H 123 (-B; T)| < ||a w A u\\ L 2 \\L Bl fix A ib 2 ^2 A ^Z^lUa 

= Ct(^J {lbiVx A iB 2/ u 2 A tB 3 /i 3 ) A *(a Bi //i A i fla // 2 A tB 3 /i 3 )) 2 

= C r (J iSxHSaHSslV)* (7.2) 
< C r E 2 (B)l, 

where E 2 (B) = f s3 \B\ 2 . To observe (1) first note that for any pair of forms: a G Q k (M), 
and /3 G Qi(N), on Riemannian manifolds M and N, on the product M x N we have 

*MxJV (ttX/(«) A tt^(/3)) = (* M 7r* M a) A (* Ar vr^/?), (7.3) 

where ttm : M x N — > M and ttn : M x N — > N are the natural projections, (the proof is 
a simple calculation in an orthogonal frame of the product and is left to the reader). Now, 
step (1) in (17. 2p follows by applying (17.31) to the integrand, and observing in the coframe 

HY- 

LBiVi A *L Bt fii = (air/ 2 A 773 - a 2 i]{ A rf 3 + a 3 i]\ A rf 2 ) A (a x ri\ - a 2 rf 2 + 03773) 
= {a\ + a 2 2 + a\) r]\ A 4 A 77* = \Bi\ 2 r]l A rf 2 A 77* = \B { \ 2 ^ 

where Bi = (ai,a 2 ,a 3 ). Step (2) in (17. 2p follows from Fubini Theorem. 

Next, we aim to provide an estimate for Cr- For this purpose we review some basic 
L 2 -theory of the operator d~ l (i.e. inverse of the exterior derivative d). The main goal is 
to estimate an L 2 -norm of the potential q u of uj in (14. 2p . Following the standard elliptic 
theory of differential forms, [31] , the potential in ( 14.21) can be obtained via a solution to 
the Neumann problem for 2-forms on T (see Appendix [Cj 

inf ~ (7.4) 

= ndcpN = 0, on dT, 
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where n stands for the normal component of a differential form, and A = d5 + 5d, 5 = ±*d*, 
(c.f. [34J ) . As T is a domain with corners we replace it by a slightly larger domain T in (S 3 ) 3 
with the same topology (i.e. T is a deformation retract of T) but with smooth boundary dT . 
(One may argue that it is not really necessary, since T is Lipschitz and elliptic problems, 
such as (17. 4j) are well posed on Lipschitz domains, [31]). Because of (i) in Theorem 14.31 we 
may use the restriction of 

« w = 5<p N , (7.5) 
to T (see Appendix Ofor justification of (17.51) ). Associated to (17.41) is the Neumann Laplacian 

A N :H 2 Q 2 N (f)^L 2 Q 2 (f), 

which has a discrete positive spectrum {X^n} and eigenvalues satisfy the variational principle 
called Rayleigh-Ritz quotient, pTJ. The first (principal) eigenvalue X\ t i$ niay be expressed as 

2 i IIJk,„l|2 



Xx,n = inf ( + v e H^ 2 N {M) n H%{M)~ 

V \m\h 

We denote the inverse of by 

G N : L 2 n 2 (f) — H 2 n 2 N (f), 

which restricts to a compact, self-adjoint operator on L 2 . As a result the spectrum of Gn 
is discrete and given as {l/Ayv}. Note that based on these considerations we may define 
d- 1 ■= 5G N . 

Theorem 7.1 (Energy bound). For every volume preserving vector field B which has an 
invariant unlinked domain T , the L 2 -energy of B on S 3 is bounded below by the third order 
helicity H 12 3(-B;T), as follows 

\H 123 (B;T)\<C T (E 2 {B))I. (7.6) 
Also, we may estimate the constant Cf. 

^<^IMIU 2( T)(MS 3 )) 3 , (7.7) 



where A^jy is the first eigenvalue of the Neumann Laplacian on Q 2 (T). 
Proof. We estimate 



\<Xu\\h = \\$4>n\\ 2 L 2 < 1 1 #w | lis + ||50jv||i2 



\(A(j) N ,(f) N )\ < \\uj\\l 2 \\(Pn\\l2, 
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where we used the Green's formula [34, p. 60] and boundary conditions of (17.41) in the second 
identity. Now, because 4>n = G^u, and it is a well known fact that ||Gtv||l 2 = XTT7' (^tv is 
compact self adjoint on L 2 ) we obtain: 



\(*u\\l 2 < nr 1MU 2 



As a result we estimate Cr'- 



□ 



Cr = \\0Cuj A L0\\ L 2 < ||^||L°°||ttu;||L 2 

1 .. .. 

— 7t IMIl 2 IM1l°°- 



Notably, the best energy estimate so far has been obtained by Freedman and He, [15J, 
for the L 3 / 2 -energy of B, in the case when B admits an invariant domain T modeled on 
an n-component link L = {L 1; . . . ,L n } in IR 3 . Their estimate is based on the asymptotic 
crossing number and reads 

1/4 ( n \ 

E 3/2 (B) > (^) ( ^ac(L fc) L)|Flux(S fc )|J ■ mmJlFlux^,)!}, (7.8) 

where the asymptotic crossing numbers ac(Lfc, L) for Borromean links can be estimated below 
by a smallest genus among surfaces in M 3 \ {L\ U...Ulj.U...U L n } with a single boundary 
component Lj. Since L 3/,2 -energy of B bounds the L 2 -energy, inequality (17. 8p leads to a 
lower estimate purely in terms of fluxes and topological data. It is not clear to the author 
if this approach can be extended to the case of invariant handlebodies considered in Section 
HI A different, more optimal estimate, has been obtained by Laurence and Stredulinsky, via 
the Massey product formula, in [26], but the proof is provided only in a special case of the 
vector field B. 

Contrary to these lower bounds, which are given in terms of topological data, the estimate 
in (17.61) depends on the geometry of the domain T, and also ||o;||loo. Unfortunately, u blows 
up on the diagonals A C (S 3 ) 3 , and as a result the estimate is meaningless when the 
handlebodies % get close to each other during the evolution of the magnetic field B. At 
this point, we need an assumption for % to stay 1cm apart during the evolution. Another 
drawback is that Ai ( jv is a geometric constant which is altered during the evolution as well. 
If we consider the situation in which the boundaries of % are invariant during the evolution, 
the estimate may be useful. Under such assumption, which occurs whenever the velocity 
field v of plasma in (II. ip is tangent to d%, the bound in (17.71) stays constant. 
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8 Comparison to the known approaches via Massey 
products. 

In several prior works [EH ES, E7j helicities were developed via the Massey product for- 
mula for /ii23- These approaches are equivalent to the one presented here in the sense that 
invariants obtained this way measure the same topological information. Most notably the 
work [25] provides an explicit expression for the third order helicity of the Borromean flux 
tubes, where the ergodic interpretation in the style of Arnold's asymptotic linking number 
is also provided. In [27] one finds the following formula for the third order helicity 



where A4 = d~ 1 (LB i fi)- This formula is valid for three distinct vector fields Bi on a closed man- 
ifold M. For invariant domains with boundary (18.11) defines an invariant provided y4j| aM = 0, 
but this only happens in certain situations (e.g. M is simply connected, and Aj's are appro- 
priately chosen). 

The most commonly known formula directly related to the Massey products was developed 
by Berger [5], in the case of Borromean flux tubes T = 1\ U T 2 U T 3 



In [5] it is expressed as a volume integral over T by applying gauge fixing. When % are 
topologically solid tori there exists a single Massey product < a±, a 2 , a 3 > in the complement 
S 3 \ T, represented by the 2-form Ai A F 23 + F 12 A A 3 . When % are handlebodies there are 
multiple Massey products, but the formula (18.21) should still be valid. So far, such extensions 
have not been considered in the literature and the volume integrals over T may be harder 
to obtain in such a case. One may also point out that ergodic interpretations of Massey 
products are more involved [25] comparing to the approach presented in Section [51 
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Appendices 

A Equations for the fro zen-in- field forms 

Given a volume preserving vector field B on M and a path t — > g(t) G SDiff(M), let 

B l := g*(t)B . 

Then for B° = B, by definition 

B t = f t B t \ t=0 =-C v B = -[V,B], (A.i; 



and as a result 



Next, calculate (dt + C-y = C 



—i B tfi\t=o — — L [V,B}^ 



dt+V, 



(d t + L v )i B tfi = dt(L B t/j) + C V L B tn 

= L B tfl + t B t(C-vf) + t'lV.B*]^ 

where in the second identity we applied the general formula: L[a,b] — C-a^b — ^b£a, and in 
the third equation the fact that V is volume preserving i.e. Cyfi = 0. As a result of (1A. 1[) 
we obtain 

(d t + C v )L B tfi = . 
Next, we justify Formula (15.51) . First use (15.61) to calculate 

L Bi (fij A fi k ) = (iBiHj) A fi k - fij A (i Bi Hk) 
l Bi (/ii A fj. 2 A /i 3 ) = (lb^i) A fi 2 A fi 3 - fix A l Bi (fi 2 A fi 3 ) (A. 2) 

= (iBi^l) A fl 2 A fi 3 - fix A i Bi fi 2 A fi 3 + fix A fi 2 A L Bi fl 3 , 
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since iB^j = for i ^ j only one term in the above expressions remains for each i. Set 
a := tB 2 t B 1 P, P ^ ^ 3 ((>S' 3 ) 3 ); since a A H\ A fi 2 A /i 3 = and a is a 1-form we obtain 

= L B:j {a A /ii A fi 2 A yU 3 ) 
= (<.s 3 Q!) A /ii A /i 2 A /i 3 - a A /ii A yU 2 A i.b 3 // 3) 

where in the last step we used ( 1A.2I) . Therefore 

{lbzLb^b-lP) A tlx A /i 2 A /i 3 = (L B2 t Bl P) A /ii A yu 2 A tB 3 // 3 . 

Analogously, {lb 2 ^b 1 P) A /ii A /12 = (iBi/3) A fi\ A ie 2 /i 2 and A H\ = {3 A iB x \i>\ which 
justifies Equation (15.51) . 

B Zero contribution of short paths to the time average 

It is clear when/3 6 fi 3 (T) is at least a C 1 on T cS 3 x5 3 xS 3 , then / = (3(B U B 2 , B 3 ) is 
continuous on T and 



/ = 

(/ +/ )(/ +/ )(/ +/ )/■ 

After expanding, it is obvious that we must show the following (for all choices), when T — > 00: 

1 

1 

1 

J>3 



/ — 0, (B.l) 



f —> 0, (B.2) 



f — ► 0. (B.3) 

a{^{x,T),x) J u{^ 2 (y,T),y) J cr(<l> 3 (z,T),z) 



Since the lengths of the short paths in S are bounded by a common constant <i, ( 1B.1I) -( |B~3~! ) 
follow immediately, e.g. for ( IB. If) we have 

/ / f\<^d(T + d)(T + d)||/|| 00 _o. 

C Notation in Section [7] 

We adopt notation from the elegant exposition in [31]. Let M be an orientable manifold 
with smooth boundary. 
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Q k (M) = C°°(M, A fc ), smooth differential forms on M. 

Q%(M) = {(f) G Q k (M) | rxf) = 0, nd<p = 0} the subspace satisfying the Neumann 
boundary conditions, (n denotes a normal component of a form along dM). 

The L 2 -inner product on Q k (M) is defined as 

(^) L 2 = / UJ A*T], 

L 2 Vt k (M), L 2 -differential forms on M. 

H k N {M) = {A G H l Vt k (M) | d\ = 5X = 0,nA = 0} the subspace of the Neumann 
harmonic fields. 

Next we justify (17.5ft . first observe that for any 7 G i/ 1 r2 fc_1 (M): 

(dj,\) L2 = [ t 7 A*nA = 0, VAG^(M), (C.l) 

where t, and n stands for respectively tangent and normal to dM components of the form. 
As a result, if uj G T-C^^M) 1 - we obtain a solution (p to the Neumann problem: 

5d(f) + d5(j) = uj, => u - dScj) = 5d(f) . (C.2) 

Formula ( 1C.1I) implies: (u — dScj)) G ^^(M)- 1 -, moreover n(u; — dScj)) = n(5d<f)) = 5n(dcj)) = 0, 
by the boundary condition in (17.4j) . If a; is a closed form, uj — d5<ft is also closed, and clearly 
coclosed by (1C.2I) . Thus uj — d5<ft is a harmonic field with zero normal component, and 
therefore it has to be in 7i k N (M), and therefore the zero form. This yields 

uj = dS(j) . 

As a result we obtain a necessary and sufficient conditions for uj to be exact: 
(i) du = 0, 

(ii) (uj, A) L 2 = 0, for all A G Tt%(M). 
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